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1. Introduction 


The Baryon Acoustic Oscillations (BAO) are a relic from the time when photons and baryons were 
coupled in the early Universe, and constitute a preferred clustering scale in the distribution of matter 
on cosmological scales. Since the physical scale of the oscillations can be inferred from observa¬ 
tions of the cosmic microwave background (CMB), one can use the measured apparent size of the 
BAO as a ‘standard ruler’ in the late-time Universe. Used as a distance measure in this way, the 
BAO are one of the most powerful cosmological observables that can be derived from large-scale 
structure (LSS) surveys - making it possible to accurately reconstruct the geometry and expansion 
history of the Universe, while being remarkably robust to systematic errors. In combination with 
the CMB and other observations, precision BAO measurements are able to decisively answer fun¬ 
damental questions about the nature of dark energy (especially its possible evolution with redshift), 
possible modifications to General Relativity, and the spatial curvature of the Universe. 

The SKA will be able to measure the 
BAO in two ways (Fig. 1). The first, us¬ 
ing the intensity mapping (IM) technique, 
seeks to reconstruct the cosmological den¬ 
sity field in three dimensions over unprece¬ 
dentedly large volumes, using the 21cm spin- 
flip transition of neutral hydrogen to trace the 
density field in redshift space. The second 
is a more conventional galaxy redshift sur¬ 
vey, which will detect many millions of dis¬ 
crete sources over a wide range of redshifts, 
with accurate redshift determinations com¬ 
ing from the 21cm line. Of the two, IM is 
the newer and less-tested method - galaxy 
surveys have been used to great effect as 
probes of LSS over the past several decades, 
whereas the first IM experiments are only 
just coming online, and are subject to a num¬ 
ber of potential foregrounds and systematics 
that are yet to be fully quantified. Neverthe¬ 
less, an IM survey is the best prospect for 
doing transformative cosmology with Phase 
1 of the SKA - surpassing all existing BAO 
constraints at low redshift, and strongly com¬ 
plementing future surveys at higher redshift. 

In this chapter, we will review the 
physics of the BAO, its suitability for use as 

a standard ruler, and the current state of the ait in BAO measurements from optical and Lyman-a 
surveys. We will then describe two possible methods for reconstructing the BAO with large surveys 
on the SKA, and present forecasts of the capabilities of SKA Phase 1 and 2 for each. 



k [Mpc 


Figure 1: Constraints on the BAO ‘wiggles’, 
/bao(^) = P{k)/Pvef{k) - 1 (where Pref(^) is a BAO- 
free reference spectrum), combined over all available 
redshift bins for three SKA galaxy/IM surveys. Full 
forecasts are given in Section 4. 
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2. BAO as a standard ruler 

In the early Universe, baryons and radiation were tightly coupled together via interactions with free 
electrons, by Coulomb attraction and Thomson scattering respectively. The tendency of overden¬ 
sities in the baryon distribution to collapse under gravity was resisted by the increase in radiation 
pressure of photons collapsing along with them, setting up acoustic oscillations in the coupled 
photon-baryon plasma. The speed of sound in the fluid gives rise to a preferred scale - the sound 
horizon, Vg - corresponding to the distance a sound wave could have travelled since the Big Bang. 
As the Universe cooled and the relative proportion of radiation to matter fell, Thomson scatter¬ 
ing became inefficient and the photons and electrons decoupled, allowing the radiation to stream 
away. With the radiation pressure gone, the sound waves stalled, and the pattern of over- and 
under-densities that they had set-up were left frozen into the baryon distribution with a typical cor¬ 
relation length corresponding to the sound horizon at this time. Cold dark matter, being coupled 
to the baryonic matter gravitationally, was drawn into the same pattern of fluctuations, leaving a 
preferred distance scale in the total matter distribution. 

One can measure this distance scale in a statistical manner by reconstructing the matter cor¬ 
relation function, i^(r) = (5M(x)5M(x-|-r)), from surveys of large-scale structure using galaxies 
or other tracers of the matter density field. The BAO scale appears as a relatively broad ‘bump’ 
feature on ~ 150 Mpc scales, and so large survey volumes are needed to measure it well. One 
only measures a relative distance scale from the correlation function; retrieving physical distances 
depends on knowing the comoving sound horizon (i.e. the physical scale of the oscillations) as 
well. This can be calibrated from the CMB angular power spectrum, for example. 

The BAO can be measured separately in the radial and tangential directions. If one makes no 
attempt to separate the two (i.e. by taking an angular average), the distance measured is d{z) = 
fsizd)/Dv{z), where Zd is the redshift of the baryon drag epoch, and the dilation scale is given by 

Dv{z)=[{\+zfDl{z)^^\ ( 2 . 1 ) 

The expansion rate, H{z), and angular diameter distance, Da{z), describe the geometry of the 
Universe, so measurements of Dy (z) over a range of redshift can be used to constrain the evolution 
of the equation of state of dark energy, w(z), for example. As shown by Eq. (2.1), however, one only 
measures a combination of H{z) and Da{z), and so model degeneracies can arise. The degeneracy 
can be broken by separating the BAO in each direction; the quadrupole of the correlation function 
is sensitive to the combination 


F{z) = {\+z)Da{z)H{z)/c, (2.2) 

and the combination of Dy and F uniquely determines Da and FI. In practise, one can separate 
out the tangential and radial directions by marginalising over the shape of the clustering pattern, 
without needing to understand the effects of redshift space distortions (Anderson et al. 2013). 

2.1 Constraining dark energy with distance measurements 

A measurement of Fl{z) and Da{z) can be directly related to the dynamics and geometry of space 
time. Specifically, the Friedmann equations relate the expansion rate to the energy content and 
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curvature of the universe via 

^ [Pm(1 + Pr(1 + Pde(^)] ~ 

where Pm and p^ are the energy density in matter and radiation today, K is the curvature of space 
and Pde(z) encapsulates all the unaccounted “dark” energy which may be contributing to the ex¬ 
pansion of the universe. This dark energy can take a wide range of different forms: scalar fields 
(or quintessence) which are analogous to the Higgs field of fhe sfandard model (Copeland el al. 
2006), gravilalional degrees of freedom lhal are a remnanf from higher dimensions, or even more 
radical modificalions lo Einstein’s General Theory of Relalivify (Clifton el al. 2012). Given ils 
exotic nalure and fhe fad lhal we cannof accounl for fhe dark energy wilh any of fhe olher, known 
forms of energy. Ibis is a unique opporlunify to explore physics beyond fhe sfandard model and 
learn somelhing qualilalively new aboul fhe fundamenlal laws of nalure. 

To learn aboul dark energy using BAO, we allempl to pin down ils time evolution. Given 
lhal we are solely focusing on Ihe dynamics and geomelry of Ihe background spacetime, Pde(^) 
is completely determined by ils value today, Pde(O), and ils equation of slate, w{z), such lhal 
d\npY)E/d\nz = 3[1 + w(z)]. Il is by characterising w{z) thal we hope to learn aboul Ihe underlying 
physics lhal is driving Ihe expansion of Ihe Universe al late times. In Iheory w{z) can have a wide 
range of behaviours, bul in practice il is useful to consider simpler paramelrizalions. A particularly 
popular and useful paramelrisalion is to - effectively - Taylor expand Ihe equation of slate in terms 
of the scale factor a = (1 +z)^^ so that 

w~ wo + Wa(l-a), (2.3) 

where wq and Wq are constants. While clearly oversimplifying the time evolution of w{z) and in 
principle only accurate for z <C 1, this parametrisation turns out to capture the behaviour of a very 
broad class of models for dark energy. Furthermore, it can be shown that, under general conditions, 
there are reasonably tight consistency relations between wq and Wa that can be used to select out 
particular types of physical model. In this paper we will use wq and Wa to quantify the precision of 
our constraints, although more general approaches have been used which involve decomposing the 
full history of w{z) into suitable functional components (Zhao et al. 2014). 

A precise measurement of wq and Wa can lead to profound insights. Current constraints are 
weak - especially from z > 1 - but are consistent with wq = — 1 and Wa = 0, which is what one might 
expect from a cosmological constant, in itself quite a remarkable result. If future constraints show 
that Wo / — 1 and w^ / 0 then it would correspond to the discovery of a new degree of freedom in 
the Universe: a new fundamental field or modifications to general relativity on cosmological scales. 
If one were to find wo < — 1 (that is, if the equation of state were to cross the “phantom divide” at 
late times (Caldwell et al. 2003)), then the consequences would be quite dramatic and one would 
have to consider a substantial revision of the current rules of field theory and gravitation. 

2.2 Systematic effects and density field reconstruction 

One of the major advantages of using the BAO to measure distances is their robustness to systematic 
effects. Since the distance information is primarily contained in the location of a feature in the 
correlation function (or equivalently the power spectrum), rather than its amplitude or detailed 
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shape, there is little dependenee on diffieult-to-ealibrate quantities like the overall normalisation of 
the matter power speetrum. Nevertheless, the non-linear growth of strueture and seale-dependent 
effeets ean blur and bias the measured loeation of the BAO feature, so must be taken into aeeount. 

Consider how galaxies respond to the growth of large-seale strueture, for example. As they fall 
toward their loeal elusters and superelusters, the large-seale galaxy pairs whieh eneode the pristine 
baryon aeoustie seale are shifted to smaller or larger separations, broadening the baryon aeoustic 
peak and deereasing the aeeuraey with whieh it ean be reeovered. The teehnique of “density-field 
reeonstruetion” (Eisenstein et al. 2006) performs an approximate eomputation of these motions 
using the observed distribution of galaxies, and henee restores objeets to the near-original position 
in the linear density field. This has fhe effeef of sharpening fhe aeousfic peak and signifieanfly 
improving disfanee measuremenfs. 

The pofenfial level of improvemenf depends on fhe number densify of fhe fraeer, as well as 
fhe individual noise of fhe realisation, buf for sample varianee-limifed surveys ean be a faefor of 
~ 2 (Seo & Eisensfein 2007). The reeonsfruefion feehnique has been applied wifh sueeess fo BAO 
measuremenfs in fhe largesf galaxy redshiff surveys af infermediafe redshiffs: fhe Sloan Digifal Sky 
Survey (SDSS) Euminous Red Galaxy sample (Padmanabhan ef al. 2012), fhe Baryon Oseillafion 
Speefroseopie Survey (BOSS) (Anderson ef al. 2013) and fhe WiggleZ Survey (Kazin ef al. 2014). 

Non-linear sfruefure formation ean also infroduee seale-dependenf effeefs. One fypieally 
fhinks of fhe presenf-day densify field beeoming signifieanfly non-linear only on seales below 
~ 30 Mpe, far below fhe BAO seale al ~ 150 Mpe. Mode eoupling on small seales, and fhe 
fael fhaf fhe observed eorrelalion funefion is a eonvolulion over all Eourier modes, means lhal non¬ 
linear effeefs have an impael on eonsiderably larger seales, however. Galaxy bias also affeels fhe 
shifls inlrodueed by non-linearifies. While non-linear effeefs lhal ehange fhe shape of fhe eluslering 
pattern ean be marginalised wilhouf biasing fhe loealion of fhe BAO fealure, fhe eombinalion of all 
effeefs ean infroduee a modesl (buf non-negligible) ~ 0.1 — 0.2% bias after reeonsfruefion (Seo el 
al. 2010;Mehla ef al. 2011). 

2.3 Previous measurements 

Sinee fhe firsl signifieanl defeelion of fhe baryon aeouslie peak in fhe SDSS Euminous Red Galaxy 
(ERG) sample by (Eisensfein el al. 2005), BAO measuremenfs have improved along wifh fhe in- 
ereasing volume mapped by galaxy redshiff surveys, and now provide per-eenl level disfanee and 
expansion measuremenfs al various eosmie epoehs. The mosl preeise eurrenl measuremenfs derive 
from BOSS (DRll, (Anderson el al. 2013)), which now offers 2% and 1% disfanee measuremenfs 
al z = 0.32 and z = 0.57 respectively, using nearly one million galaxies covering 8,500 deg^ of sky 
and a volume of 13 Gpc^. BOSS has also measured BAOs in fhe sfruefure of fhe Eyman-a foresl on 
fhe sighllines lo quasars, providing a 2% disfanee measuremenl al z = 2.34 (Delubac el al. 2014). 

Measuremenfs al lower precision have previously been reported by fhe WiggleZ Dark Energy 
Survey (~ 4% measuremenfs in Iwo independenf bins al z = 0.44 and z = 0.73, (Kazin el al. 2014)), 
fhe 6-degree Eield Galaxy Survey (4.5% measuremenl in fhe local Universe al z = 0.1, effeclively 
serving as an independenf determination of Hq (Beutler et al. 2011)), and the final SDSS ERG 
sample (Padmanabhan el al. 2012). BAOs have also been delecled in pholomelric redshiff surveys 
al intermediate redshiffs (Seo el al. 2012). 
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Upcoming optical and near-infrared redshift surveys enabling BAO measurements include the 
extended BOSS project (eBOSS), which will provide precision distance measurements in the range 
z < 1, the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), targeting 2 < z < 4, the 
Euclid and WEIRST satellites, and proposed ground-based projects such as the Dark Energy Survey 
Instrument (DESI) and the VISTA/4MOST telescope. Ref. (Eont-Ribera et al. 2013) describe how 
a combination of these facilities should provide percent-level distance constraints over the range 
0 < z < 4. In this chapter we discuss whether SKA surveys are able to compete in this landscape. 

3. Experiment and survey design 

There are two main methods that can be used to measure the BAO with the SKA: a galaxy redshift 
survey and an intensity mapping survey. In this section, we describe the relative merits of the 
two methods, and define baseline specifications for both types of survey on the Phase 1 and 2 
configurations. 

3.1 HI galaxy redshift survey 

Galaxies are biased tracers of the cosmological density field. By detecting many individual galaxies 
and measuring their redshifts, one can constrain the matter correlation function (or equivalently, the 
power spectrum). As discussed above, redshift surveys at optical wavelengths have been used to 
great effect for cosmology, and with the SKA one will be able to do the same in the radio. 

Most important is the choice of target galaxy population. In particular, one requires galaxies 
with an easily measurable emission/absorption feature for measuring redshifts - the most appropri¬ 
ate for the SKA being HI (Santos et al. 2014b). The high spectral resolution of SKA receiver sys¬ 
tems rivals the redshift resolution of optical spectroscopic surveys, and their large bandwidths allow 
a wide redshift range to be covered in principle. The Phase 1 arrays have insufficient sensitivity to 
yield competitively-sized samples of HI galaxies, however - even with optimistic assumptions, a 
10,000 hour survey over 5,000 deg^ with SKAl-MID or SUR will achieve an RMS flux sensitiv¬ 
ity of Srms ~ 70 — 100 /iJy, equating to roughly 5 x 10^ galaxies out to z ~ 0.5 for a 5a detection 
threshold. This is worse than the expected yield from the full 10,000 deg^ BOSS survey, which will 
be completed long before Phase 1 sees first light. SKA2, on the other hand, will be far more sensi¬ 
tive, reaching Si-ms ~ 5 /iJy for a 10,000 hour survey over 30,000 deg^, even with a more stringent 
10a threshold. The expected yield for such a survey is ~ 10^ galaxies between 0.18 < z < 1.84, 
far surpassing any planned optical or near-infrared survey for z < 1.4. The predicted number den¬ 
sity and bias of HI galaxies for SKAl-MID (including MeerKAT dishes), SKA 1-SUR (including 
ASKAP dishes), and SKA2 are given in Table 1; more details on the sensitivity calculation are 
given in Santos et al. (2014b). 

There are a number of potential systematic effects that can affect galaxy redshift surveys. In 
large optical galaxy surveys at least, stars are a major contaminant - while one can distinguish stars 
from galaxies by their colour, bright stars effectively mask galaxies behind them, leading to a com¬ 
plicated angular selection function on the sky. This is less of an issue in the radio, although other 
contaminants, such as diffuse galactic synchr'otron emission and non-galaxy point sources, can also 
cause problems for the source-finding algorithms used to compile the galaxy catalogue. Source- 
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SKA 1-MID 


Zc 

n{z) [Mpc 

b{z) 

0.05 

2.92 X 10-^ 

0.678 

0.15 

6.74 X 10^3 

0.727 

0.25 

1.71 X 10^3 

0.802 

0.35 

4.64 X 10-4 

0.886 

0.45 

1.36 X 10-4 

0.975 


SKAl-SUR 

Zc 

n{z) [Mpc-^] 

b{z) 

0.05 

4.14 X 10-2 

0.664 

0.15 

8.00 X 10-2 

0.724 

0.25 

1.56 X 10-2 

0.802 

0.35 

3.39 X 10-4 

0.890 

0.45 

7.86 X 10-2 

0.989 

0.55 

1.91 X 10-2 

1.099 

0.65 

4.77 X 10-^ 

1.221 

0.75 

1.23 X 10-^ 

1.357 

0.85 

3.21 X 10-2 

1.507 


SKA2 

Zc 

n{z) [Mpc-2] 

b{z) 

0.23 

4.43 X 10-2 

0.713 

0.33 

2.73 X 10-2 

0.772 

0.43 

1.65 X 10-2 

0.837 

0.53 

9.89 X 10-2 

0.907 

0.63 

5.88 X 10-2 

0.983 

0.73 

3.48 X 10-2 

1.066 

0.83 

2.05 X 10-2 

1.156 

0.93 

1.21 X 10-2 

1.254 

1.03 

7.06 X 10-4 

1.360 

1.13 

4.11 X 10-4 

1.475 

1.23 

2.39 X 10-4 

1.600 

1.33 

1.39 X 10-4 

1.735 

1.43 

7.99 X 10-2 

1.882 

1.53 

4.60 X 10-2 

2.041 

1.63 

2.64 X 10-2 

2.214 

1.73 

1.51 X 10-2 

2.402 

1.81 

9.66 X 10-^ 

2.566 


Table 1: Predicted number density and bias of HI galaxies as a function of redshift, for SKA 1-MID (in¬ 
cluding MeerKAT, Band 2), SUR (including ASKAP), and SKA2 (Az = 0.1 bins). All assume 10,000 hour 
surveys, over 5,000 deg^ (SKAl-MID and SUR) and 30,000 deg^ (SKA2). The detection thresholds are 
chosen to be 5(7, 5(7, and 10(7 respectively. See Santos et al. (2014b) for more details. 


finding in radio data is also made more challenging by needing to search through 3D datacubes 
resolved in angle and frequency, rather than just sets of discrete 2D images. 

One should also be careful of source evolution effects. For example, the luminosity function 
of the tracer population is expected to change with redshift, modifying the number of galaxies that 
can be detected. Depending on the tracer, this limits the useful redshift range of a survey, and 
complicates its selection function. This evolution is commonly parametrised as an evolution of the 
galaxy bias, b{z), and the principal effect on the BAO is to change the shot noise by changing the 
effective galaxy number density, n{z). It is also possible that the bias can become scale-dependent, 
however, which introduces the possibility of systematically biasing the BAO scale. 

As discussed in Section 2.2, redshift space distortions and effects on non-linear scales also 
affect the measured power spectrum. A simple model for the impact of RSDs and non-linearities 
on the galaxy power spectrum is given by (Kasier 1987; Seo & Eisenstein 2007), 

Ptot(k) = (3.1) 

aNL(z,At) = aNLD(z) {\+f{z)^i^[2 + nz)]y (3.2) 

where D{z) is the growth factor, f{z) = dlogD/dloga is the linear growth rate of structure, and 
P{k,z) is the isotropic matter power spectrum. The leading term of Eq. (3.1) is the anisotropy due 


7 


















Measuring BAO with future SKA surveys 


Philip Bull 


to RSDs (/I = cos 6), and the exponential term models the “washing-out” of redshift information 
due to incoherent non-linear velocities, characterised by the velocity dispersion, aNL- These effects 
are discussed in more detail in Raccanelli et al. (2014). 


3.2 Intensity mapping survey 

If one is only interested in measuring the matter distribution on large scales, there is not strictly any 
need to resolve individual galaxies. Instead, one can perform a survey with relatively low angular 
resolution that detects only the integrated intensity from many unresolved sources. This is called 
intensity mapping (IM), and allows extremely large volumes to be surveyed very efficiently. 

The SKA will be capable of performing large IM surveys over 0 < z < 3 using the redshifted 
neutral hydrogen 21cm emission line (Santos et al. 2014a). Neutral hydrogen is ubiquitous in 
the late universe, residing principally in dense regions inside galaxies that are shielded from the 
ionising UV background, and the 21cm line is narrow and relatively unaffected by absorption or 
contamination by other lines. It is thus an excellent tracer of the matter density field in redshift 
space. The background brightness temperature of the HI emission is (Bull et al. 2014) 


Tt 


3 hc^Aio 
3271 kBmpV^^ 


H{z) 


^m{z)Pc,o, 


(3.3) 


where Aio is the Einstein coefficient for emission, and flHi(z) is the fraction of the critical density 
today, pcfi, in neutral hydrogen. Assuming that HI is a linearly-biased tracer of the total matter 
density field, we can relate fluctuations in the brightness temperature to the (Fourier-transformed) 
redshift-space matter density perturbation. 


5n = TMk) = + (3.4) 


where the anisotropic terms are due to RSDs and non-linear growth (Eq. 3.1). One can then 
measure the 3D redshift-space matter power spectrum, (5M(k)5^(k')) = (27r)^5(^)(k — k')P(k), 
by mapping out the brightness temperature distribution in angle and frequency, V = Vhi/(1 +z). 



SKA 1-MID 

SKAl-SUR 


Band 1 

Band 2 

Band 1 

Band 2 

Tinst [K] 

28 

20 

50 

30 

^min 

0.35 

0.00 

0.58 

0.00 

^max 

3.05 

0.49 

3.05 

1.18 

Pmin [MHz] 

350 

950 

350 

650 

Pmax [MHz] 

1050 

1760 

900 

1670 

Ddish [m] 

15 

15 

12 

12 

5v [kHz] 

50 

50 

50 

50 

Osur [10^ deg^] 

25 

25 

25 

25 

^dish ^ ^beam 

254 X 1 

254 X 1 

60x36 

60x36 


Table 2: Baseline IM survey specifications for the Phase I SKA configurations. SUR is equipped with 
PAFs; the assumed scaling of the PAF FOV with frequency is explained in Santos et al. (2014b). 
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Example survey specifications for Phase 1 of the SKA are shown in Table 2. One of the key 
decisions in designing an IM survey for the SKA will be whether to use the array in autocorrela¬ 
tion (single-dish) or cross-correlation (interferometer) mode. As one is interested in mapping out 
extended structure on comparatively large angular scales rather than detecting individual galaxies 
that subtend only small angles, the angular sensitivity is an important factor. Roughly speaking, 
single-dish experiments are sensitive to angular scales between the field of view (FOV) of a single 
dish and the full area of the survey, whereas interferometers can only see between the scales corre¬ 
sponding to their minimum and maximum baseline lengths. Since the minimum baseline can be no 
smaller than the diameter of a single dish, the maximum possible angular scale that interferometers 
are sensitive to is set by the single-dish FOV (i.e. the primary beam). The two modes are therefore 
in some sense complementary (Fig. 2). 

For the ~ 15m dishes of the SKA, 
the physical scales corresponding to the 
BAO are best matched to an autocorre¬ 
lation survey for z < 1, and an interfer¬ 
ometric survey at higher z. Having the 
BAO features fall somewhere within an 
instrument’s resolution window is only 
a minimum requirement, however - ide¬ 
ally, one should maximise its sensitiv¬ 
ity at the relevant angular scales too. 

This is trivial in the single-dish case, 
which has an essentially flat response 
over its full range of angular scales, but 
for interferometers there is a strong de¬ 
pendence on the detailed baseline dis¬ 
tribution. For IM with the SKA, a 
high density of short baselines is op¬ 
timal, suggesting an array configura¬ 
tion with highly-clustered stations. For 
the currently-proposed SKAl configu¬ 
rations, however, the density of short 
baselines is not high enough, and auto¬ 
correlation mode always wins out on sensitivity at low redshift. An interferometric survey would 
be better suited to BAO detections at high redshift, z > 1-5, but this is a less interesting range for 
constraining dark energy. A proposed dense aperture array component of SKA2 would likely be 
better suited to intensity mapping; aperture arrays have the large FOV and high density of short 
baselines necessary for sensitivity to the angular scales of the BAO at lower z. 

Though promising, the intensity mapping methodology is not yet mature, and it remains to 
be seen whether a number of potentially serious issues can be overcome. The most obvious of 
these is the presence of foreground emission from the galaxy and extragalactic point sources, 
which strongly contaminates the HI signal. Fortunately, despite being several orders of magnitude 
brighter, most of the foregrounds are spectrally smooth and so can be readily distinguished from all 



k j_ [Mpc ] 

Figure 2: Maximum and minimum transverse (angular) 
scales that can be probed as a function of redshift in single¬ 
dish (autocorrelation) and interferometer modes, for a rep¬ 
resentative SKA array (15m dishes with baselines restricted 
to a 1 km core, and 25,000 deg^ survey area). The BAO 
(dashed vertical lines) and comoving horizon scale (grey 
band) are shown for comparison. (N.B. the array will also 
be sensitive to structure in the radial (frequency) direction.) 
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but the largest-scale modes of the fluctuating cosmological signal. More insidious is the leakage of 
polarised foreground emission into the total intensity channel; this is not spectrally smooth because 
of frequency-dependent Faraday rotation caused by the partially-incoherent magnetic field of the 
galaxy (Alonso et al. 2014). The issue of foreground contamination is explored in more detail in 
Wolz et al. (2014). 

Another significant issue for single-dish surveys in particular is correlated (“1//”) noise. A 
well-known problem in CMB analysis, correlated noise from the receiver system and elsewhere 
dominates the uncorrelated white noise component on long timescales. This causes the signal-to- 
noise ratio to grow more slowly as a function of integration time, and introduces striping artefacts 
into maps of the emission, seriously hindering the recovery of large-scale modes of the signal. For 
the CMB, this is handled by using receivers with low knee frequencies (< 1 Hz), scanning rapidly 
across the sky, and filtering out timescales longer than l//knee from the time-ordered data. An 
SKA autocorrelation survey would have to adopt similar strategies to make a wide-angle survey 
viable. Finally, ground pickup/spillover is also an issue for autocorrelation surveys, and must be 
suppressed (e.g. with ground shielding), or mapped out and subtracted. 

4. Forecasts for the SKA 

In this section we present forecasts for BAO distance measurements - and the resulting constraints 
on cosmological parameters - for several SKA configurations, and compare them with what will 
be possible with other methods in around the same timeframe. The forecasts are based on the 
Fisher forecasting formalism developed in (Bull et al. 2014), using the number counts from Table 
1 (galaxy survey), and experimental specifications set out in Table 2 (IM survey). Complementary 
forecasts for redshift space distortions are given in Raccanelli et al. (2014). 

While there are a number of possible combinations of arrays, frequency bands, and sur¬ 
vey modes for both Phase 1 and 2, for compactness we have chosen to concentrate on the best¬ 
performing configurations. For example, for Phase 1 we neglect Band 1 of SUR and MID for 
galaxy surveys, as they will detect comparatively few HI galaxies, and we consider only autocor¬ 
relation mode for intensity mapping, since the angular resolution of an interferometric IM survey 
on the SKA is poorly matched to the BAO for z < 1. We present forecasts only for a galaxy survey 
for SKA2, although it should be noted that an IM survey on a mid-frequency aperture array should 
be able to provide similar constraints on the BAO out to z Ri 2. 

4.1 Fisher forecasting 

To accurately characterise the expected performance of a given experiment, one would ideally per¬ 
form a full simulation, incorporating a variety of potential systematic and instrumental effects, and 
running the simulated data through the actual analysis pipeline. This is computationally intensive, 
and detailed aspects of the hardware and signal may not yet be known, as is the case here. Fisher 
forecasting takes a simpler approach, instead using the expected properties of the signal and noise 
for an experiment to derive a Gaussian approximation to the likelihood for a set of parameters to 
be measured. Though clearly not definitive, Fisher forecasts at least take into account important 
effects like correlations between parameters, and are sufficiently accurate for understanding the 
relative performance of different experiments. 
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The key procedure in Fisher forecasting is to construct the Fisher matrix, F, for a set of pa¬ 
rameters, 6. When inverted, this yields an estimate of the covariance matrix for the Gaussianised 
likelihood. The Fisher matrices for both IM and galaxy redshift surveys can be written in the form 


1 f d^k , , 


d 'T d T 

j^logC .j^logC 


(4.1) 


where is the total covariance of the measured signal, consisting of the true sig¬ 

nal (S) and noise (N), and Veff(k) = Vphys (C^/C^)^ is the effective volume of the survey (which 
covers a physical volume Vphys)- Only is a function of the cosmological parameters of in¬ 
terest, {6}. For a galaxy survey, the total signal is the galaxy power spectrum plus shot noise, 
= Ptot(k) + \/n{z). For intensity mapping, the expression is more complicated: oc 

and oc rgyg5'area/hot®(k), where B is a window function that depends on the angular and fre¬ 
quency resolution of the radio telescope (full expressions for both interferometer and autocorre¬ 
lation experiments are given in Bull et al. (2014)). We have neglected the effects of foreground 
subtraction here. 

To calculate the Fisher matrix, we must specify a fiducial cosmology and a set of parameters to 
be measured (including nuisance parameters). We adopt the Planck best-fit ACDM model (Planck 
Collaboration 2013), 


h = 0.67, Ha = 0.684, flK = 0, Q.b = 0.049,>v = -1, n, = 0.962, ag = 0.834, N^n = 3-046, 

(4.2) 

and in the first instance forecast for the parameters 

{a_L =DA(z)|fid./-DA(z),a|| =H{z)/H{z)\M.,(y%f{z),o%b{z),ns,Om.}, (4.3) 

where ‘fid.’ denofes evaluation in fhe fixed fiducial cosmology. Since we are only inferesfed in 
fhe BAO here, we discard information on fhe disfance paramefers {a^, a||} from all ofher sources. 
In practise, fhis means fhaf we only keep fhe ferms wifh derivatives of fhe BAO pari of fhe power 
speclrum wifh respecl fo a, i.e. df'QAoik)/da, where we have splif fhe isolropic power specfrum 
info a smoolh pari and an oscillating pari, P{k) = [1 -|-/BAo(^)]7’smooth(^)- We also assume fhaf no 
reconslrucfion of fhe densily field is performed. 

Finally, one can projecl from fhis sel of paramefers info anolher fhaf corresponds direclly lo 
fhe parameters of a cosmological model. We will consider fhe projeclion of fhe dislance/growlh 
functions, {a^, a||,/}, info fhe parameters 

{/l,n^,OA,O^,W0,Wa,aNL,C^8,^(z)}- (4.4) 

This corresponds lo filling a cosmological model lo fhe disfance measuremenfs exfracfed from fhe 
BAO, and fhe growfh rate measured from fhe anisolropy of fhe correlafion funclion. 

4.2 BAO forecasts for the SKA 

Predicted constraints on the BAO feature in the power spectrum for three different SKA surveys 
were shown in Fig. 1, as a function of scale. For Phase 1, it is clear that an intensity mapping 
survey will have much better overall sensitivity to the BAO than a galaxy survey when constraints 
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Figure 3: BAO-only fractional constraints (68% CL) on the expansion rate (top panel) and angular diameter 
distance (bottom panel) for SKA IM and galaxy redshift surveys. The bias has been marginalised as a free 
parameter in each redshift bin. The degradation of the Da constraints at higher z for the IM experiments is 
caused by their falling angular resolution (c.f. Fig. 2). 


are combined over the full redshift range. A galaxy survey with SKA2 will far surpass both of 
these (subject to various systematic effects). 

Fig. 3 shows forecast constraints on the expansion rate and angular diameter distance for the 
proposed SKA surveys, in bins of width Az = 0.1 (galaxy surveys) or Av = 60 MHz (IM surveys). 
Forecasts for a Euclid galaxy survey are also shown for comparison, based on the predicted number 
counts and bias in Amendola et al. (2013). For Phase 1, a galaxy survey will not be competitive with 
other BAO measurements owing to insufficient sensitivity. An IM autocorrelation survey will be 
significantly more powerful, providing constraints on H{z) that are similar to the (sample variance- 
limited) Euclid experiment, but over a significantly wider redshift range - Band 2 of both SUR and 
MID will yield sub-2/3% constraints beyond z — 2. Constraints on the angular diameter distance 
will be considerably worse at higher redshift due to the limited angular resolution in autocorrelation 
mode, however (c.f. Eig. 2). A galaxy survey with SKA2 will be sample variance-limited over 
30,000 deg^ for 0.4 < z < 1.3, surpassing all other planned surveys over that range. ^ 

* An IM survey on a Phase 2 dense mid-frequency aperture array operating from 450 MHz upwards could provide 
similarly tight constraints on Da and H all the way out to z = 2 if a large enough collecting area could be built. 
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Figure 4; Forecast constraints on {wo,Wa) for several SKA configurations and Euclid, in combination with 
Planck and BOSS. All other parameters have been marginalised, including CIk, and the bias is free per z bin. 


Expected constraints on the dark energy equation of state parameters (Eq. 2.3) are shown in 
Eig. 4 for a subset of these configurations, combined with forecast Planck CMB and BOSS low 
redshift BAO constraints.^ Eorecasts for a future Euclid galaxy survey are shown for comparison, 
also with Planck+BOSS included. While unable to match Euclid on raw figure of merif^ (EOM = 
69, versus 129 for Euclid), fhe Phase 1 IM survey has a highly complemenfary redshifl range and 
wide survey area (wifh pofenfially ~ 100% overlap), and should be complefed in fhe early 2020’s, 
making if fhe low redshifl dalasef of choice for joinl analyses wifh Euclid and olher high redshifl 
experimenls like WEIRST and DESI. Combining surveys will be imporlanl for pinning-down fhe 
nalure of dark energy in fhe medium lerm, as “Slage IV” galaxy surveys may nol offer sufficienl 
precision on Iheir own lo discriminate belween mosl dark energy models (Marsh el al. 2014). In 
Ihe longer term, Ihe galaxy survey wilh SKA2 will be able lo achieve a subslanlially larger EOM 
of around 310. 

5. Conclusions 

The Baryon Acoustic Oscillations imprinl a dislance scale into Ihe dislribulion of mailer on large 
scales lhal can be used to conslrain Ihe expansion history and geomelry of Ihe Universe. In torn, 
Ihis can be used to conslrain key cosmological parameters, potentially allowing us to answer fun- 
damenlal questions aboul Ihe nalure of dark energy such as whelher il evolves wilh time. The BAO 

^We have assumed that the BOSS sample is statistically independent from the surveys that we combine it with. 

^The dark energy figure of merit is defined as FOM = l/.^det(T^' Iwo.wo) (Albrecht et al. 2006). 
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are remarkably robust to systematie effeets, and have been sueeessfully measured to high preeision 
by a number of optieal galaxy redshift surveys and Lyman-a forest observations. 

The SKA will be able to measure the BAO through two different types of survey: HI galaxy 
redshift surveys, where millions of individual galaxies are deteeted and their redshifts measured 
from the HI emission line; and HI intensity mapping surveys, where the integrated HI emission 
from many unresolved galaxies is used to reeonstruet the eosmologieal density field on large scales. 
Redshift surveys are a tried and tested technique, but require high sensitivity to capture enough 
galaxies. Intensity mapping has yet to mature, but is potentially a much more efficient way of 
detecting the BAO. Both suffer from a number of potential systematics, such as those associated 
with the non-linear evolution of the cosmological density field on small scales. 

Phase 1 of fhe SKA will be able fo produce compefifive consfrainfs on fhe (mosfly radial) BAO 
af redshiffs relevanf for dark energy if a large IM aufocorrelafion survey can be performed. The 
Phase 1 arrays lack fhe sensifivify fo defecf enough galaxies fo produce inferesfing consfrainfs from 
a redshiff survey, and fhere are an insufficienf number of shorf baselines fo make an inferferomelric 
IM survey worfhwhile (unless higher redshiffs are fargefed). Nofe fhaf an IM aufocorrelafion survey 
on an early deploymenf (pre-Phase 1) sysfem may also be able fo produce useful consfrainfs on fhe 
expansion rafe so long as sufficienl survey lime (~ 10,000 hours) can be oblained. 

SKA2, on fhe olher hand, will have fhe sensifivify fo produce an immense galaxy redshiff 
survey over almosl 3/4 of fhe sky, surpassing all olher planned BAO measuremenls al 0.4 < z < 
1.3. This should allow if fo pin down fhe equalion of slale of dark energy wilh unprecedenled 
precision, assuming various syslemalic effecls can be overcome. A mid-frequency aperlure array 
could achieve similar precision al higher redshiff (oul fo z ~ 2) if a sufficienlly large and dense 
collecling area can be used, bul we have nol considered Ihis possibilify in delail here. 

Finally, while we have concenlrafed on fhe BAO as fhe mosl robusl dislance measure, redshiff 
space dislorlions and even fhe overall shape of fhe power speclrum conlain a greal deal of exlra 
information lhal can be used fo consfrain dark energy. In Ihis sense, fhe forecasls here represenl 
fhe mosl conservalive eslimales of fhe cosmological consfrainfs lhal can be achieved wilh fhe SKA. 
If sufficienl conlrol over syslemalics can be achieved, considerably lighler measuremenls of w(z) 
and fhe growlh of slruclure can be expecled by using Ihese olher probes loo (Raccanelli el al. 2014). 
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